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We study an atomic Bose-Fermi mixture with unpolarized fermions in an optical lattice. We 
obtain the Mott ground states of such a system in the hmit of deep optical lattice and discuss the 
effect of quantum fluctuations on these states. We also study the superfluid-insulator transitions of 
bosons and metal-insulator transition of fermions in such a mixture within a slave-rotor mean-field 
approximation, and obtain the corresponding phase diagram. We discuss experimental implications 
of our results. 
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I. INTRODUCTION 

Recent experiments on ultracold trapped atomic gases 
have opened a new window onto the phases of quan- 
tum matter^'. A gas of bosonic atoms in an optical or 
magnetic trap has been reversibly tuned between super- 
fluid (SF) and insulating ground states by varying the 
strength of a periodic potential produced by standing 
optical wavesi*^. This transition has been explained on 
the basis of the Bose-Hubbard model with on-site re- 
pulsive interactions and hopping between nearest neigh- 
boring sites of the lattice i^i^i^i^ii. Further, theoreti- 
cal studies of bosonic atoms with spin and/or pseu- 
dospin have also been undertaken^i^ i^^'^^ . These studies 
have revealed a variety of interesting Mott phases and 
superfluid-insulating transitions in these systems. On the 
fermionic side, the experimental studies have mainly con- 
centrated on the observation of paired superfluid state&H 
and the BCS-BEC crossover in such systems near a Fes- 
hbach resonanceiS,. 

More recently, it has been possible to generate mix- 
tures of fermionic and bosonic atoms in a trapiiii^. Ini- 
tially, the main focus of such experimental studies were 
to generate quantum degenerate Fermi gases, through 
sympathetic cooling with bosons. However, a host of 
theoretical studies followed soon, which established such 
Bose-Fermi mixtures to be interesting physical systems 
in their own righti^iiiiiSiii^i^S, exhibiting exciting Mott 
phases in the presence of an optical lattice. In all of 
these works, the spin of the fermions in the mixture is 
taken to be frozen out due to the presence of the mag- 
netic trap. However, more recent works reported in Refs. 
[21I and [12, considered a Bose-Fermi mixture in an optical 
trap, where the spins of the fermions can be dynamical 
degrees of freedom^^ . It has been shown in Ref. '2T that 
the interaction between the bosons and the fermions in 
such a mixture can enhance the s-wave pairing instability 
of the fermions. 



In this work, we consider a Bose-Fermi mixture in an 
optical trap and in the presence of an optical lattice 
and study the Mott phases and the metal/superfluid- 
insulator transition for the fermions/bosons of such a 
mixture using a slave-rotor mean-field theory^!. The mo- 
tivation for such a study is two-fold. First, it has been 
shown in Ref. [isl that the Mott phases of Bose-Fermi mix- 
tures in a magnetic trap are interesting in their own right. 
In the present study, we chart out the Mott phases of the 
Bose-Fermi mixture in an optical trap, where the Fermion 
spins are dynamical degrees of freedom, in the deep lat- 
tice limit for a wide range of parameters. As expected, 
we find that the corresponding Mott phases obtained are 
much richer than their counterparts studied in Ref. ITsl . 
Second, the metal/superfluid- insulator transition of the 
fermions/bosons in such an interacting Bose-Fermi mix- 
ture has not been studied before. Here we develop a self- 
consistent slave-rotor mean-field theory to study such a 
transition for the Mott phases which do not have density- 
wave order with broken translational symmetry and use 
it to obtain at least a qualitative understanding of the 
effect of interaction between the fermions and the bosons 
on the metal/superfluid- insulator transition. 



In what follows, we shall assume that the atoms are 
confined using an optical trap so that Fermions spins are 
not frozen out. We shall, however, ignore the effect of 
the harmonic trap potential which is a standard approxi- 
mation used extensively in the literature^^. The starting 
point of our study is the Bose-Fermi Hubbard Hamilto- 
nian that has been developed earlier, with similar ap- 
proximation regarding the trap potential, from underly- 
ing microscopic dynamics of the atoms in the presence of 
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Here c^o- is the fermionic destruction operator with spin 
cr at site i, bi represents bosonic destruction op- 

erator at site i, n^^^^ denotes Fermion(Boson) number 
operators, ii?(B) and ijlf{b) are nearest neighbor hopping 
matrix elements and chemical potentials for the fermions 
(bosons), Ubb and Uff are the on-site Hubbard re- 
pulsion for bosons and fermions respectively, and U fb 
denotes the relative interaction strength between the 
bosons and the fermions. In what follows, we shall take 
the bosons and the fermions to have fixed chemical poten- 
tials Hb{f) and same on-site repulsion Ubb = Uff = U 
and consider A = Ufb/U and ij = tF/tB as parameters 
which can be freely varied. The justification of this choice 
is briefly outlined in Sec. [Ill Further, we shall only deal 
with case of a square bipartite lattice in this work since 
this is simplest to realize experimentally. 

The organization of the rest of the paper is as fol- 
lows. In the next section, we identify the Mott phases 
of Eq. [TJ Next, in Sec. Illli we introduce the slave rotor 
formalism and use it within a mean-field approximation 
to study the metal/superfluid-insulator transition of the 
Bose-Fermi Hubbard model (Eq. [T|) . This is followed by a 
discussion of possible experiments in Sec. IIVI A compar- 
ison of the Mott-Hubbard phase diagram obtained using 
the projection operator technique with those obtained 
from mean-field theories:^'^^'^ and standard strong cou- 
pling expansions^ - is given in Appendix [Xj 



II. MOTT PHASES 

In this section, we chart out the Mott phases of the 
Bose-Fermi system. To do this, we first obtain the phases 
of the system in the Mott limit {Ib = ^f — 0) and 
then obtain fluctuation corrections over these states to 
0(i|(^)/{/2). 

Before obtaining the Mott phases for the Bose-Fermi 
mixture, let us look briefly into the parameters of the 
Hubbard model (Eq. [1]) . These can be determined from 
the microscopic quantities such as the potential depths 
Vf{b) due to the laser seen by the atoms and their recoil 



energies Ep^^^ 



t^k\/2raF(B) where is the wave- 
vector of the laser and mF{B) are the masses of the 
fermions (bosons). The potential depth seen by the atoms 
depend on the detuning of the laser from their natural 
wavelengths Xf{b) of the fermions (bosons). In fact, it 
can be shown that, the ratio of the lattice potentials seen 
by the fermions and bosons are^i 
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where AXf{b) — Xl — Xf{b) denote the detunings for the 
fermions (bosons) and ^f{b) are the corresponding natu- 
ral linewidths. Since the ratio of the natural linewidths is 
generally close to unityi^i^, we see that one can tunc the 
ratio of the lattice depths seen by bosons and fermions 
by varying the detunings. 

In terms of these quantities, we have^i 
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where apFi olbBi and apB are the s-wave scattering 
lengths for interaction between two fermions, two bosons 
and a Fermion and a Boson respectively. These scatter- 
ing lengths also can be varied either by choosing different 
species of fermions or bosons or by tuning them using 
Feshbach resonance. Further, as we have discussed be- 
fore, by choosing the laser detuning we can also make the 
fermions and bosons see either similar or very different 
lattice potentials. Therefore, instead of calculating these 
parameters from the microscopies, we shall aim to por- 
tray a general picture of the Mott phase diagram. Since 
the experimental possibilities are limitless, for the sake of 
brevity, we choose Ubb ~ Uff — U and vary the ratios 
A — Ubf/U and rj ^ tF/tB- It is clear from the above 
discussions that such a situation can be always achieved 
in experiments. We shall consider some such specific ex- 
amples in Sec. IIVI 

Next, we consider the Bose-Fermi Hubbard Hamilto- 
nian in the Mott limit. In this limit, the on-site states 
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where we have scaled all energies by U and = 
f^F{B)/U- It can be seen from Eqs. [71 that two states 
\nQ,nQ — l) and — l,n^) are degenerate when 

l)(n^ - n^) ^ f^F~ fJ-B 



(A 



whereas three states j? 
IjTIq > are degenerate when 
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The conditions of these degeneracies, of course, depend 
on our choice of parameters of the model. It is also to be 
noted that in Eqs.[H]and[ni tiq and tiq in the ground state 
are themselves functions of fj,'p, fj,'^ and A, and have to be 
determined by minimizing Eq.[7]subject to the constraint 
of 71^ and being integers. 

The ground state phase in the Mott limit {tf, = tp = 0) 
diagram can be obtained by numerically minimizing the 
ground state energy (Eq. [7]) for integers < < 2 
and . For the sake of brevity, we carry out the nu- 
merical computation for /i'^ = /i^ = /i and present the 
phase diagram as a function of /i and A. The phase di- 
agram for A > is shown in Fig. [T] We find, as ex- 
pected from the results of Ref.[21), for A > 1, the fermions 
and the bosons repel each other out from a given site so 
that a site is occupied by either a Boson or a Fermion, 
but not both. Such states were dubbed as "compos- 
ite" states of bosons/fermions with a correlation hole of 
fermions/bosons in Ref . '2T. However, in the present sce- 
nario, the spins of the fermions are dynamical degrees of 
freedom which leads to richer variety of possible phases, 
as we discuss below. 

As expected from discussions leading to Eq. [H 
a large portion of phase diagram has degenerate 



ground states corresponding to 



1q = 0, = 1) and 
|n^ = l,n^ = 0). Note that in terms of the origi- 
nal fermions, the degeneracy is actually threefold corre- 
sponding to states |0, t), |0, J,) and |1, 0). This degeneracy 
is lifted by quantum fluctuations due to the presence of 
small but finite tp and is- This can lead to three differ- 
ent ground states as sketched in Fig. [2| A) an antiferro- 
magnetic state of fermions with no bosons, B) A state of 
one boson per site and no fermions, and C) A state with 
fermions and bosons being the nearest neighbors with an 
antiferromagnetic order for the fermions. The energies 
of these states can be estimated using a straightforward 
^i^%{F)/^^) perturbation theory and are given by 



Ea 



Ef 



Ec 



Nztj, _ 
if " 
2Nzt% 

Nzt% 
' 2XU 



Nzt\ 
U ' 



(1 + ^') 



(10) 



where N is the total number of sites in the system and 
z is the coordination number of each site. Comparing 




FIG. 1: (Color online) Ground state phase diagram in the 
atomic limit for /i — fi'g — ^'p and A > 0. The phases are 
marked by values of (n^,no') in the ground state. For large 
A, the system tries to avoid putting bosons and fermions at 
the same site. The threefold degeneracy between (1,0) and 
(0, 1) as well as (2, 1) and (1, 2) occur for a large portion of 
the phase diagram. In addition, there are doubly degenerate 
regions such as (2, 0) and (0, 2). These degeneracies are lifted 
by virtual hopping process for small but finite tg and tp 



the energies of the states from Eq. [TOl we find that 
state A is favored over B and C when rf' > 2 and 
A > (1 + ?7^^)/2. Similarly the state B is favored for 
rf <2 and A>(l+7y^)/4 and state C for small A when 
(1 -I- rf )/{2\) > Max [2, 77^]. The corresponding phase 
diagram is shown in Fig. [2l 

With only a nearest-hopping term in the Hamiltonian, 
the antiferromagnetic ordering in the ground state C is 
in fact frustrated. Indeed, Ath-order (in t) virtual hop- 
ping leads to an antiferromagnetic coupling both between 
next-nearest-neighbor sites and next-next-nearest sites 
(Fig. [3]). In practice, however, this frustration is sup- 
pressed by 0(i'^/J7) corrections that arise from the next- 
nearest neighbor hopping t' (not included in the Hamilto- 
nian ([2])), leading to an antiferromagnetic order at wave- 
vector (tt/o, 0) (or (7r/\/2a, 7r/-\/2a) in a reference frame 
tilted by 45°), in two dimensions as shown in Fig. [21 
This should be contrasted with the usual (7r/a,7r/o) or- 
dering realized either for state A or for the (1,1) state 
shown in Fig. [H Additional hopping amplitudes (e.g. 
the next-next-nearest-neighbor hopping t") are expected 
to be smaller and will not affect the antiferromagnetic 
ordering of the ground state C. A similar consideration 
applies for the degenerate states (2, 1) and (1, 2) where 
quantum fluctuations will similarly lift the degeneracies. 
Notice that the possibilities of having a ferromagnetic 
state of fermions where all sites are uniformly occupied 



FIG. 2: (Color online) Possible ground states A, B and C 
that results from lifting the classical degeneracy of the ground 
states {riQ ^Uq) = (1,0) and (0,1) due to quantum fluctua- 
tions. Notice that the antiferromagnetic ordering of state C 
is different from state A and is a consequence of very weak 
next nearest neighbor hopping. Similar states will result when 
the degeneracy of the states (2, 1) and (1, 2) is lifted. 



FIG. 4: (Color online) Ground state phase diagram in the 
Mott limit for A < 0. The number of fermions per site is 
Uq — 2 while the number of bosons per site is shown in each 
phase for large negative A. 



FIG. 3: Virtual hopping (dashed lines) leading to frustra- 
tion of the antiferromagnetic ground state C. This frustration 
is lifted and an antiferromagnetic ground state stabilized by 
virtual hopping generated by the kinetic coupling t' between 
next-nearest neighbor sites (located at opposite corners of the 
square lattice unit cell). 



by T or J, fermions never occur since such a state neces- 
sarily suppresses Fermion hopping due to PauU principle 
and is thus higher in energy compared to the state A. 

Finally, we point out that lifting of degeneracy by 
quantum fluctuations also occurs for the twofold degen- 
erate states labeled as (2, 0) or (0, 2). Note that the state 
(2, 0) corresponding to two fermions per site requires that 
all second order virtual hopping processes are suppressed, 
unless higher bands are involved. Hence such a state 
is energetically unfavorable and is never realized. The 
two other states are D) a homogeneous state with two 
bosons per site and E) a state with alternate arrange- 
ments of two bosons and two fermions per site. The 
energies of these states are given by Ed — —6Nzt^/U 
and Ee = -Nzt%{l + rf)/{2X - l)U . The latter state 



(E) is thus favored over state (D) for if > 6 (2 A — 1) — 1. 

In contrast for A < {i.e. Ubf < 0), there are no 
degeneracies. The (1, 1) state persists at small negative 
A. For |A| > 1, the system always has two fermions per 
site with uq bosons, where hq is the integer which mini- 
mizes i5[no] — ~M"o + no{nQ — l)/2 — 2|A|no. The phase 
diagram is shown in Fig. [4] as a function of fi and A. The 
phase diagram corresponds to Mott phase of bosons cou- 
pled with non dynamical fixed number of fermions. Note 
that analogous composite states with uq bosons and one 
Fermion were found and dubbed as " composite Fermion" 
states in Ref. [111. Here we have almost identical states 
at large negative A with the difference that there are two 
fermions per site (instead of one per site as in Ref. [2ll ) 
owing to the fact that Fermion spins are not frozen in 
the present study. 



III. SLAVE ROTOR MEAN FIELD THEORY 



In this section, we construct a slave-rotor mean-field 
theory for the coupled Bose-Fermi problem and use it 
to study the metal/superfluid- insulator transition in this 
system. We develop the formalism for the mean-field 
theory in Sec. IIII Al and discuss the results obtained in 

secnn] 



5 



A. Formalism 

We begin by implementing the slave rotor formalism^^^ 
in the present context. This key observation behind this 
formalism is that the fermionic Hubbard Hamiltonian 
(Eq. [2]) can be mapped onto a Hamiltonian of free auxil- 
iary Fermions coupled self-consistently to a quantum ro- 
tor. The chief advantage of this representation is that the 
quartic interaction term of the original fermionic Hub- 
bard model can now be represented by a quadratic term 
in the rotor variables and can thus be treated exactly in 
the Mott limit. This feature makes this technique suit- 
able for studying Hubbard models in the strong-coupling 
regime. Further, as shown in Ref. [23, the metal-insulator 
transition of the original Fermions can be looked upon as 
the order-disorder transition of the rotors which facili- 
tates the study of metal-insulator transition. 

To begin with, we note the identity 
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F F 1 \ ^ F 
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so that, up to a constant term, one can write Eq. [2]as 
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Next, following Ref. |2J, we introduce the slave rotor rep- 
resentation for the fermions. The key observation here 
is that the spectrum of Hp in the Mott limit depends 
only on the total Fermion number rua which can be 
represented by eigenvalues of angular momenta of a 0(2) 
rotor. Thus we write the physical fermion annihilation 
operator as 



(13) 



where fia denotes the annihilation operator for the 
pseudo-fermion and 9i denotes the rotor variable at site 
i. The corresponding angular momentum of the rotor is 
denoted by Li = —idg-. In this representation, a physical 
fermion state can be written as a product of the pseudo- 
fermion state and a rotor state as 



|ci,C2...Cq)^ = |/l,/2, •■•/q)J^ = Q - 1) 



(14) 



Here I denotes the eigenvalues of the rotor angular mo- 
mentum L and |ci,C2...cq) is the antisymmetric combi- 
nation of Fermionic states for Q = I + 1 fermions at site 
i. Note that whereas for ordinary rotors I can take all 
possible integer values, here it is constrained to range 
between —1 and 1 by the operator identity 



L, 



(15) 



With the following constraint, one can write Eqs. [T2land 
m in terms of the rotor variables as 



Uf = ^tp J2 {flfj- exp [i {9. - 9j)] + h.c) 



M^-|)E/U. + |E^^' (16) 
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with Li related to J2afLf^<y by Eq. [TH In the Mott 
limit (ti? = tB = 0), we can implement the constraint 
(Eq. [T5)) exactly and this procedure leads to Eq. [7] with 
fiQ = Q = ^ -|- 1 in a straightforward manner. Note also 
that the quartic interaction term J/X^i^^T'^ii Eq. [2] 
has now been replaced by a quadratic term U J2i ^1/2 in 
Eq. [161 This has been done at the expense of generating 
a non- linear coupling between the auxiliary Fermions and 
the rotors, as is evident from the first term of Eq. [TBI 

When tp,tB 7^ 0, the above mentioned constraint con- 
dition can not be implemented exactly and we need to re- 
sort to mean-field approximation. The slave-rotor mean- 
field theory for the present system can be developed by 
a straightforward generalization of the formalism devel- 
oped in Ref. [H. To begin with, we write the system 
Hamiltonian Ti. as 



n 

Hf ^ - 



Hr + Hf + Hb (18) 
- 5^ COS {9, - 9,) + ^L^ + hL,{19) 

<ij> i 
<ij>,a 

+ E(-A'F + U/2 + h + \UnB)flha (20) 
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{fib - XUnp)'Y'' 



(21) 



where we have implemented the constraint Eq. [TSj us- 
ing an auxiliary field hi which has been replaced by 
it's saddle point value h at the mean- field level. Here 
we also treat the coupling between the bosons and 
the fermions within mean-field approximation by re- 
placing the Fermion/bosons density operators nf^^ by 
their averages ub/f- This amounts to replacing the 
term Y.i>^Unfnf by Y^i^^nffiB in Eq. ^ and by 
Xj XUnffip in Eq. [21] Within this mean- field approx- 
imation, the coupling term acts as a density-dependent 
shift in the chemical potentials for the bosons and the 
fermions. The effective hopping matrix elements tf^ and 
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rfl^ in Eqs. fT9l and [20l are given by 



(22) 



H, 



with the assumption that averages such as (exp(0i — 6j)) 

and (j2a fiafjc^ ^''"S '^SS'I ea,ch bond^i. 

The next step is to approximate the rotor model by an 
effective single site model 



i ^ 

K = -2^ 



U . 
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(23) 
(24) 



Note that this approximation makes sense only when the 
Mott ground state does not have density wave order of 
any kind either for bosons or fermions. Using the fact 
that under this approximation the becomes a single 
site Hamiltonian, we define 



Z = {cos{9))l^ 



(25) 



and use Hf to compute all average involving the 
Fermionic fields. This yields 

K = 4 (cos(6')>^^ y" deenpiZe- fiF + h + XUub) 
= 4(cos(0))„ e (26) 



(L) = 2 

= 2nF - 1 



J deD{e)e {Ze - ^iF + U/2 + h + XUhb) 



- 1 
(27) 



where the density of states I?(e) is defined as D{e) = 
J d'^kSie- ek)/i2Tr f, and Cfc = ~2tJ2^^^.^y,,cos{ki) is 
the kinetic energy of the fermions (all momenta are mea- 
sured in units of lattice spacing), and e and uf are the 
average fermionic kinetic energy and density respectively. 
Comparing the expression of uf with that for the free 
fermions uf = fl'l^ deD{e), where is the chemical po- 
tential for the free fermions at T = 0, one has the relation 

Z^o = fiF -U/2 - h - XUfiB (28) 
(L) = 2nF-l (29) 

Notice that since Z vanishes at the transition at ^f = 
IJ,p, one gets fip — U/2 = h + XUfiB- Thus Z acts as the 
order parameter for the metal-insulator transition of the 
fermions. 

Eqs. HSl m HZl US] and [H have to be self-consistently 
solved to obtain the ground state of the system. However, 
to do this, one needs to obtain the ground state of the 
bosonic Hamiltonian (Eq. I2ip and compute the average 
value of the boson density Tlb- Since at the mean-field 
level, the average Fermionic density np enters the boson 



Hamiltonian as a shift in the chemical potential, one can 
use the projection operator technique developed in Ref. 
[m for obtaining ub- It was shown in the context of 
two-species Bosons that the projection operator method 
compares well with quantum Monte Carlo results^^. 

The projection operators for the boson Hamiltonian 
fEg. [21]) can be constructed following the procedure of 
Ref. [ll|. It is given by 

Pi = {\n^)^{nEl)®[\n^)^{nE\) (30) 

where is the Boson occupation per site for the Mott 
ground state and I denotes the link connecting two neigh- 
boring sites i and j. 

The hopping term for the bosons can be rewritten in 
terms of sum over links as 



T 



b]h. 



h.c. 



(31) 



where i and j are near neighbor sites connecting the link 
I. In this notation, one can now divide the hopping terms 
into two parts 



Ti 



{PiTi + TiPi) + Pi^TiPi^ (32) 



where P^ = 1 — P; . It is then easy to see that the term 
= Tl acting on the ground state takes one out of 
the ground state manifold. The idea is therefore to seek 
a canonical transformation operator S which eliminates 
Ti to 0{tB/U) from the low energy effective Hamilto- 
nian: [iS^Hq] = —T^, where Hq denotes all the on-site 
terms in Eq. The effective low energy Hamiltonian can 
be obtained by the usual Schrieffer- Wolff transformation 
method 



(33) 



1 



= ifo + T° + T] + - [tS, [tS, Ho]] + ... (34) 

Note that this is equivalent to a systematic tB/U expan- 
sion and all the omitted terms denoted by ellipsis are at 
least 0{t%lU^). 

The next task is to find out the canonical transforma- 
tion operator S in terms of the projection operators Pi. 
Following Ref. [ill, we guess the form of the S to be 



S = iaY,[PuTi] 



(35) 



where the coefficient a is to be determined by the con- 
dition [iS^Hq] = — T^. To do this we use the operator 
identities 

[Pm,Ho] - UPiTi, [TiPi,Ho\^~UTiPi (36) 
and evaluate [iS, Hq\ to be 

\iS,Ho] - 



a 



J2{PiTi+TiPi) = -T'- (37) 
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Thus we find that setting a = U we obtain the expression 
for 5* = i J21 [P^Ti] /U, which ehminates Ti to 0(^5 /C/) 
from the low energy effective Hamiltonian. 

The effective Hamiltonian can be now rewritten by 
substituting the condition [iS, Hq] = —T^ in the last term 
of Eq.lMlas 

H* = Ho + T°-J2 [[PhTi] , + Ti]/2] (38) 
1,1' 

With some algebra we now reach the final form of the 
low energy effective boson Hamiltonian which takes into 
account all t'^/U^ fluctuations 



H* 



Ho + r 



1 

u 



J2 {PiT^Pi - TiPiTi] 
I 

J2 \ PiTiTv - TiPiTi' 
1,1' I 



PiTiTuPi, ~ TiPiPi,Ti, 



(39) 



where / and I' are nearest neighbor links. One can now 
use a on-site variational wavefunction in the same way as 
in Ref. [O 

= l[{a\nl^),^ + b\ni + l)^ + c\ni -1)^) m 

i 

to minimize the ground state energy Eq — {^v \ H* Y^v) 
and obtain the corresponding boson density fiB and su- 
perfluid order parameter A 

fiB = |aG|'7i^ + |6G|'K + l) + |cG|'K-l)(41) 

A = ^n^ + la*abG + \f^c*(jaG (42) 

where aQ, be and cq are values of the coefficients a, 
b and c in the variational ground state. A compari- 
son of the phase diagram obtained by minimizing Ey = 
H* I'iiy) for A = with analogous phase diagrams 
obtained from mean-field theory^i^ and defect phase cal- 
culations to 0{t\/U'^)^ is presented in App. \K\ 

Eqs. gT] and combined with Eqs. [IHl [H HZl US] 
and can now be solved self consistently to obtain 
a mean-field description of the coupled Fermi-Bose sys- 
tem near the metal/superfluid- insulator transition points 
which are signaled by the onset of a non-zero Z or 
A. This method, therefore, allows for a self-consistent 
treatment for the coupled Bose-Fermi problem near 
the metal/superfluid- insulator transition point, provided 
that the insulating ground state preserves translational 
symmetry. 



B. Results 

In this section, we discuss the results of application of 
the formalism developed in the last section to the prob- 
lem at hand. Here we shall concentrate on the (1, 1) 
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FIG. 5: (Color online) Plot of the order parameter A and Z 
for n/U = 0.7 and 77 = 5. Non-zero values of A / Z signals 
superfluid/metal-insulator transition for the Boson/fermions. 
The symbols are as follows; black squares and green solid line 
(A and Z respectively for A = 0), red circles and cyan dotted 
line (A and Z respectively for A = 0.3), and blue triangles 
and magenta dashed line (A and Z respectively for A = 0.5). 
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FIG. 6: (Color online) Plot of the average densities ub and 
np for the bosons and Fermions for ^/U = 0.7 and rj = 5. The 
deviation of the densities from their quantized values in the 
Mott state signals the onset of metal-insulator transition for 
fermions and superfluid-insulator transition for the bosons. 
The symbols are as follows: black square and red circle (fis 
and Hf respectively for A = 0), green uptriangle and blue 
downtriangle (ns and np respectively for A = 0.3), magenta 
dotted line and cyan solid line (n_B and np respectively for 
A = 0.5). 
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FIG. 7; (Color online) Plot of A and Z for ^/(7 = 0.4, 
?7 = 4 and A = and A = 0.3. All symbols have the same 
meaning as Fig. (5] The plot serves as an illustration that re- 
pulsive interaction between bosons and fermions can enhance 
metal/superfluid transitions. 

Mott state for the following reasons. First, since we 
would be interested in studying metal-insulator transi- 
tions for fermions together with the superfluid-insulator 
transition of bosons, we would like to concentrate on 
Mott states which has one Fermion per site. Second, 
in the formalism developed in Sec. IIII Al we treat the 
Bose-Fermi interaction term within mean-field approxi- 
mation, we would like to restrict ourselves to the param- 
eter regime |A| = \UfbIU\ < 1, where the mean-field 
results are expected to be more accurate. 

To demonstrate the effect of interspecies interaction, 
we first concentrate on a fixed value of = Ms — 
0.7/7 and rj = tp/tB — 5, and study the onset of 
metal/superfluid- insulator transition as a function of 
tp/U for a few representative values of A. A plot 
of Z and A for this case, is shown in Fig. [5] while 
the fermionic and bosonic densities are plotted in Fig. 
[SI The results of these plots can be understood as 
following. Consider fixing rj and gradually increas- 
ing tp so that the fermions/bosons moves towards a 
metal/superfluid- insulator transition point. As long as 
the bosons/fermions are in the Mott state, their densi- 
ties are pinned to Uq = 1 or = 1, and hence the 
fermions/bosons see a fixed chemical potential f^^pf^g^ = 

/i — XUtiq (tiq). For our chosen value of rj, the metal in- 
sulator transition occurs before the superfluid-insulator 
transition of the bosons at tp{X) — tp {fi°p)- Notice that 
tp{X) is a non-monotonic function of A for a given /i and 
77. Hence the metal-insulator transition for the fermions 
can either be enhanced (A = 0.5) or hindered (A = 0.1) 
due to the Bose-Fermi interaction. Once the fermions 
have delocalized, the density of fermions changes with tp 



for a fixed chemical potential n as seen in Fig. [6l Hence 
the effective chemical potential seen by the bosons now 
becomes a function of both ts = tp/rj and A. Thus by 
increasing ts, we actually traverse a curve with a finite 
slope in the ii—tB plane in contrast to the non-interacting 
(A = 0) case. Consequently, the superfluid-insulator 
transition of the bosons occur at i^(A) which can be quite 
different from ^^(A = 0). We note this effect may lead 
to both enhancement [^^(A) > t^(A = 0)] or hindrance 
[ig(A) < ^sCA = 0)] of the superfluid-insulator transition 
of the bosons depending on the chosen values of rj and /i, 
as can be seen by comparing Figs. [5] and [Tj For the choice 
of At/;7 = 0.7 and 77 = 5, we find that t%{\) > t%{X = 0), 
whereas the reverse case is realized for n/U = 0.4 and 
rj = A. In the latter case, both the metal-insulator tran- 
sition for the fermions and the superfiuid-insulator tran- 
sition for the bosons are enhanced for A > 0. Hence 
we conclude that a repulsive Bose-Fermi interaction can 
either enhance or hinder the onset of metal/superfluid- 
insulator transition of a coupled Bose-Fermi mixture held 
at a fixed chemical potential. Notice that this effect is 
absent if the boson and the fermion densities are held 
constant individually, since in that case, there is no in- 
fiuence of the fermionic metal-insulator transition on the 
bosonic superfiuid-insulator transition within the mean- 
field theory. This is of course an artifact of the present 
mean-field approximation. Clearly, the dynamics of the 
fermions/bosons should play an important role in the 
transition. For example, near the metal-insulator of the 
fermions, there will be density fiuctuations which give the 
bosons a chance to hop to a neighboring site even if they 
would remain localized in the absence of such fiuctua- 
tions. A treatment of this effect requires analysis of the 
slave-rotor model beyond the present mean-field approx- 
imation and is outside the scope of the present study. 

Finally, we present a plot of the critical hopping 
strength in Fig. [5] as a function of rj and A for a 
fixed representative ji/U = 0.7. We note that at small 77, 
the superfiuid-insulator transition takes place when the 
fermions are still in the Mott state with their density 
pinned at rPp — 1. Consequently, the transition for the 
boson, within the simple mean-field theory, is the same as 
that occurring for /ics — fJ- ~ At/; becomes maximum 
when /ieff/C/ ~ 0.4 or A « 0.3. However when 77 is large, 
the fermions have already undergone the metal-insulator 
transition when the boson are at their transition point 
and hence have (np) ^ 1. Consequently the bosons see 
a different fXctt ^ fJ- — AC/ (np) and hence the value of 
changes. This is reflected in bending of the phase bound- 
ary in the right half of Fig. [H] Analogous plots for tp will 
have qualitatively same features. 

Before closing this section we would like to make a few 
qualitative comments. First, the analysis of the phase 
diagram for negative A with |A| < 1 can be carried out 
in a similar manner and one obtains qualitatively simi- 
lar results in this case. Second, at large and negative A, 
the Mott states correspond to two fermions and no (A) 
bosons localized per site. In this case, with decreasing 
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IV. EXPERIMENTS 




FIG. 8: (Color online) Plot of tg/t^ ^is a function of ?7 and A 
for n/U = 0.7. For small r), the fermions remain in the Mott 
state when the superfluid-insulator transition of the bosons 
takes place whereas for large r; they have already undergone 
a metal-insulator transition. 



lattice potential, the bosons would undergo a superfluid- 
insulator transition with the fermions still in the Mott 
phase. Therefore the superfluid-insulator transition of 
the bosons can be described, within mean-field theory, 
by a Hubbard model of bosons with effective chemical 
potential /ioff — IJ- ~ 2XU. As we decrease the lattice 
depth, the fermions would eventually delocalize when the 
hopping coefficient tp becomes comparable to the energy 
gap between the single particle levels {En — 5Eii) in 
the potential wells. However, this is a large energy scale 
{En ^ U, XU) for deep lattices and so for any reasonable 
value of 1], one expects the superfluid-insulator transition 
of bosons to occur while the fermions are still the Mott 
state. Third, we would like to note that the slave-rotor 
mean-field theory worked out here can not be applied for 
states with broken translational symmetries in terms of 
boson and fermion numbers ( such as the state C shown 
in Fig. [2]) since we have used a single-site approximation 
for the slave-rotor mean-field theory. Also the present 
slave-rotor treatment is expected to be more inaccurate 
for states with X ^ 1 since the fermion-boson interac- 
tion is treated at a mean-field level within our scheme. 
Nevertheless we note that qualitatively we would expect 
a Mott-superfluid/metal transition from these phases as 
A is reduced keeping the density of bosons and fermions 
constant. Finally, the present theory can not address the 
question of possible superconducting instability of the 
fermions in the metallic state. This issue is discussed in 
Ref. [2l|. A generalization of our mean-field description 
of the superfiuid/metal- insulator transitions to address 
such issues and inclusion of the effect of quantum fiuctu- 
ations beyond the mean-field treatment used here remain 
open problems to be addressed in future works. 



A large part of the Mott phases and the superfiuid- 
insulator transition of the Bose- Fermi mixtures discussed 
here can be experimentally accessed using standard ex- 
periments on specific Bose-Fermi mixtures. Before going 
into details of specific systems, let us first outline some 
typical experiments that can be performed on these sys- 
tems. One such experiment that is routinely carried out 
in atomic systems is measurement of momentum distri- 
bution of atoms in the trapi. This is typically done in a 
time of flight measurement by letting the atoms fly out 
by dropping both the lattice and the trap and then mea- 
suring the position distribution of the expanding atom 
cloud. Such a position distribution measurement yields 
information about momentum distribution of the atoms 
within the trap. These experiments obtain qualitatively 
different signatures for atoms in the Mott and the super- 
fluid states for bosona^*^, but can not distinguish between 
the Bose and the Fermi atoms. To achieve this distinc- 
tion, one needs to pass the expanding atom cloud through 
a Stcrn-Gerlach apparatus. Since the fermions and the 
bosons have different spins (or total angular momenta), 
they will get separated during this process and can thus 
be distinguished. Such Stern-Gerlach experiments have 
been carried out with bosonic atoms in Refs. [13 and[28l 
and their generalization to present systems should be 
straightforward. Finally, the antiferromagnetic order of 
fermions in the Mott phases (either the (1,1) phase in 
Fig. [T] or the phases A and C shown in Fig. [2]) can be de- 
termined by measuring spatial noise correlations of the 
expanding cloud in a time of flight measuremen t^^i'^'^ . 

Let us now consider some specific Bose-Fermi mixtures 
that have been realized experimentally. One such system 
is ^Li and ^Li mixtures with app = ubb — 5ao and 
ubf > '^^'^^ , where oq is the Bohr radius. The value of 
aBF has not been unambiguously measured in this sys- 
tem, but is expected to be positive indicating a repulsive 
interaction^yi. For this mixture, mp ^ niB and depend- 
ing on the value oi Ubf and by varying the frequencies 
of the laser providing the optical lattice (as discussed 
in Sec. [ni, we may realize different points on the phase 
diagram shown in Fig. [TJ Of particular interest is the 
case where qbf > app^^- In this case, for one atom per 
site, we expect to obtain the one of the Mott phases A, 
B or C shown in Fig. [2] depending on the specific value 
of A ~ ubf/o-ff and 77 = tB/tp which can be varied by 
slightly changing the ratio Vp/Vb (Eq.[Sp^. In this case, 
one can scan a large part of the phase diagram shown in 
Fig. [21 The possible antiferromagnetic orders of fermions 
in phases A or C (Fig. [2]) can be determined by measuring 
the spatial noise correlation measurements^i. The states 
B and C can also be distinguished by passing the ex- 
panding clouds through a Stern Gerlach apparatus^L^^. 
Other possible Mott phases, as shown in Fig. (TJ are also 
possible if there are more than one atom per site. These 
phases can be detected analogously. 

On the other hand, if aBF < off, one may access 
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the (1,1) phase (Fig. [T]) with one fermion and one bo- 
son per site. Here in the Mott state, we would find two 
separate clouds for bosons and fermions in the Stern- 
Gerlach measurement with an antiferromagnetic state for 
the fermions which can be deduced in the spatial noise 
correlation measurement. Further, one can now access 
the superfiuid/metal- insulator transition in this system 
by reducing the depth of the laser producing the opti- 
cal lattice. The superfluid-insulator transition can be 
directly accessed by measuring the measurement momen- 
tum distribution of the bosons. This provides us a direct 
measurement of tg. One can now also change rj as dis- 
cussed in the last paragraph and access (77) for a given 
A. This would provide access to a line of constant A in 
the phase diagram of Fig. [51 

Another Fermi-Bose mixture which has been realized 
so far ^°K-^^Rb mixture"^"*. Here one expects apF — 
aggM, and negative apB indicating an attractive inter- 
action between bosons and fermions. The magnitude of 
apB is also measured in Ref . Is^ and is found to be around 
3.6aFp, although with a large (about 40%) uncertainty. 
In this system, we expect to find Mott phases where two 
fermionic atoms sits on the same site with no{X) bosonic 
atoms. Such states can also be detected in experiments 
by passing the expanding clouds through a Stern Gerlach 
apparatus as discussed earlier. The superfluid-insulator 
transition of the bosons can also be accessed by lowering 
the lattice depth. 

In conclusion, we have studied a Bose-Fermi mixture 
in an optical lattice trapped by an optical trap. We 
have sketched a generic phase diagram for the possi- 
ble Mott states of these systems and also studied the 
superfluid/metal- insulator transition using a slave-rotor 
mean-field theory. We have also discussed definite ex- 
periments that can be performed on specific experimen- 
tally realized systems that can probe at least part of the 
above-mentioned phase diagrams. 

PM's visit to the Cavendish Lab is supported by the 
EPSRC (UK) and by a Visiting Scholars Grant from 
Trinity College. 



APPENDIX A: COMPARISON OF PHASE 
DIAGRAMS 

In this section, we compare the phase diagram for a 
single species Mott-Hubbard system obtained using the 
projection operator with those obtained using mean-field 
theories^i^i^ and strong-coupling expansions'^. To do 
this, we consider the case of zero coupling between the 
bosons and fermions (A — 0). Our starting point is the 
Bose-Hubbard Hamiltonian of Eq. [21] with A = 0. By 
tracing the same set of steps, as in Sec. IIII A[ we then 
obtain the effective Hamiltonian H* (Eq. [35]). The next 
step is to obtain the variational energy using the wave- 
function ^„ (Eq. [40]) . For the purpose of variational en- 
ergy computations, it is sufficient to consider ^E"^ with 
real coefficients a,b and c (Eg. [^D]) . This amounts to set- 



ting the phase of the superfluid order parameter A (Eq. 
[42]) to zero and does not affect the variational energy. A 
straightforward calculation then yields 



Ey 

Eq 

El 
E2 



(^y\H*\^y) ^Eo+Ei+E2 
6Epb^ + SEhC^ 

—ztBO^ [(rio + 1)^^ + n-oc^] 

_zt|no(no + l) (^4_2^2^2) 

z{z - l)t|no(no + 1) / 4 



(Al) 
(A2) 
(A3) 



U 



2z(z-l)t|Vno + l 
U 

'{2no + l)b^c^] 



[2bc {b'^ino + l) + c^no) 
(A4) 



where SEp — (— /i + C/?io) and SEh = (ju — [/(tiq — 1)) 
are the on-site energy costs of adding a particle and a 
hole respectively to the Mott phase, and z — 2d denotes 
the coordination number for a d dimensional hypercubic 
lattice. The phase diagram can now be obtained by min- 
imizing the variational energy Ey for given (Ib /U, ^/U). 
The Mott-superfluid phase boundary then corresponds 
to the minimum value of tdfJ-) for which the superfluid 
order parameter A (Eq. [35]) is non-zero. Notice that ig- 
noring the 0{t\/U^) terms amount to setting £^2 = 0. 
Our numerical results in this section, however, retains all 
terms in E2. 

Next, we obtain the expression of t™^ using mean- 
field theory. This can be done in a standard manner as 
shown in Refs. ,3n5ii6ii37i and the mean-field critical hop- 
ping strength in d dimensions can be obtained^ 



^^0 + 1 _^ 

Uno — fi ji — U (no — 1) 



(A5) 



where riQ denotes the boson occupation number of the 
Mott phase from which the phase boundary is ap- 
proached. 

Finally, we compare the phase diagram obtained from 
minimizing the variational energy E,, with the strong- 
coupling expansion developed in Ref. l26l . The main idea 
behind the strong-coupling expansion is that at the phase 
transition point, for a given /i/J7, the defect state, which 
corresponds to an additional particle or hole added to the 
Mott state, becomes energetically more favorable. Thus, 
the energy difference of the particle or hole defect states 
with the Mott state given to 0{t%/U^) hy^ 



zt^ 

= 6Ep ~ ztsino + 1) + ^ (S'^o + 4) no 



^ ''B 

u 



no {no + 1) 



zti 



SEh - ztsno + ^ (5^0 + 1) ("o + 1) 



^ ''B 

U 



no (no + 1) 



(A6) 



(A7) 
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vanishes at = ^?oriJ. The phase boundary is ob- 
tained by finding the critical hopping — Mm[t^,t'^]. 

Before resorting to numerical evaluation of the phase 
diagram from all different techniques, we would like 
to clarify the following points. First, although both 
the defect state calculations of Ref. IH and the 
projection operator technique outlined here captures 
some contributions of t'^/U'^ fluctuations, they are 
not identical to each other. To see this, we con- 
sider a second order virtual process for a defect state 
with one additional particle (hole) |no + (— )l)j \no)j 
\n + (-)2)^ \n - \n+ (-)2)^ \n)^, where i and 

j are nearest neighbor sites on the square lattice. 
This process, which after summing over all sites, gives 
0(i^/C/^) energy contributions Ep = — z^t^no("-o + 
2)/2U for particles and Eh = -zH\{nl-l)/2U for holes. 
Note that all the states involved in such a process lie out- 
side the low energy manifold and hence are not captured 
within the projection operator technique even if states 
with |no ± 2 > are incorporated in the variational wave- 
function (Eq. [40|) . On the other hand, the projections 
operator technique together with the variational wave- 
functions leads to terms in E2 (Eq. IA4|) which involves 
product of states with one additional particle and hole 
(terms which involve product of the coefficients b and 
c). These terms, which becomes important mostly near 
the tip of the Mott lobe, are necessarily absent in de- 
fect state calculations in Ref. [H which considers energy 
lowering due to a single particle or hole added over the 
Mott state. Thus we find that the best way to compare 
different approaches is to compare the phase diagrams 
obtained using them. 

Second, if we neglect the 0{t\/U'^) terms, the saddle 
point equations for h and c can be easily obtained by 
minimizing the variational energy in Eqs. IA2 land IA3l 



(no + 1) - — ^ 



no 



Ej^ 

zts 



b = 2{no + l)b^ + noc^b 

c = uqc^ + {no + l)6^c (A8) 



where we have used the constraint + b'^ + ~ 1. For 
the Mott phase, the solution to Eq. lASl is b = = c and 
a = 1. Note that this ensures that the density in the Mott 
state is pinned to (n) = no- Solutions with non-zero b 
and c occurs when in the SF phase for which ts > tc- We 
find that in the superfluid phase near the Mott transition 
line, Eq. IA8I admits the following solutions. For /i/C/ < 
no/{2no + 1) one gets tc = t'^ — SE^/zno and 



V2\ ts 



6 = 0, 



(A9) 



whereas for fi/U > no/(2no + 1) one has tc = t^ 
SEp/z{no + 1) and 



V2\ tB 



1/2 



c = 0. (AlO) 



Eqs. IA9I and lAlOl shows that near the phase transition 
b,c <ti a, which is crucial to our analysis. Further, as nor- 
mally expected in a second order quantum phase tran- 
sition, the coefficients b and c are continuous across the 
transition. Next, we discuss inclusion of 0{t\/U^) terms. 
If these terms (Eq. lA4p are included, obtaining analytical 
solutions for b and c becomes difficult as it amounts to 
solving two coupled cubic equations. However, we have 
checked during numerical evaluation of the phase dia- 
gram minimizing Eq. lAll ( which retains second order 
terms in Eq. IA4p that b and c are always small compared 
to a near the phase transition and that near the ends of 
the Mott lobe their numerical values are well reproduced 
by Eqs.[A9]and|Ai0l 

Finally, we comment about our choice of variational 
wavefunction in which we have only retained states with 
one additional particle and hole per site over the parent 
Mott state. In principle, one can retain states with two 
or more particle per sites which leads to a more compli- 
cated trial wave function. For example, one can retain 
the states |no ± 2) so that the variational wavefunctions 
become 



IV'') 



IV''),; = a|no) +fc|"o + 1) +c|?io - 1) 

+rf|no + 2) +e|no-2) (All) 



which leads to a variational energy to 0{tB/U) (Eq. 
E', = (V^'l HoiX = 0) + TO 1^') = E^, + E{^, where, 



SEpb^ 



SEhC^ 



iz/iC- + SE2pd- + SE2he 
6E2p = -2/i + C/(2no + 1) SE2h ^ 2^i - 3C/(no - 1) 



-Ztb 



{no + l)a^b^ + noa?c? 



{no 



+ \/no{no - l)c^ae + y^(n^r+"2)(n^H~T)6^arf 
+2^/{no - l){no + 2)bcde + {no + 2.)(fb^ 

+abc (d^/no{no + 2) + e^Jnl-l] (A12) 



To analyze the phase diagram obtained from Eqs. lAllI 
and IA12i we first note that in the Mott phase a = 1 
and b — c = d = e = 0, as discussed before. Next let 
us discuss the minimization of E'^ near the Mott-SF 
transition line. When d — e — 0, the MI-SF transition 
occurs at ztc = Mm{t^,tc), with 6, c ^ 1 and a ~ 1. So 
to find out whether non-zero d and/or e is favorable for 
energy minimization, we need to find the effect of turning 
on a non-zero d and/or e near this transition point and 
check whether it leads to lowering of the variational 
energy. We find that turning on a non-zero d and/or e 
near this point, as can be seen from the expression of 
E[.^j (Eq. IA12[) , leads to energy gain <C ztc (since in the 
expression of E[^, d and e always appear as a product 
with b and c which are small), whereas the energy cost 
fEq. lXTT|) is 0{SE2p) and/or 0{SE2h)- It turns out that 
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FIG. 9: (Color online) Phase diagram for d = 2 and no = 1 
calculated by minimizing variational energy Ev (black dots), 
mean-field theory (red solid line) and defect energy calcu- 
lation to 0{t%/U^) (blue solid line). The phase boundaries 
computed by minimizing i?„ and by defect energy calculations 
compare well to each other but differs from the mean-field 
theory near the tip of the Mott lobe. 
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FIG. 10: (Color online) Phase diagram for d = 3 and no = 
1. All notations are same in Fig. |9] As in 2D, the phase 
boundaries computed by minimizing and by defect energy 
calculations compare well to each other but differs from the 
mean- field theory near the tip of the Mott lobe. 
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the latter is always greater than the former which leads 
to a net non-zero energy cost. As a specific example to 
demonstrate this point, let us consider ~ C/ so that 
for tB tc = tP, b ^ a,nd c = (Eq. Then 
turning on a non-zero d results in an energy change SE — 

d {8E2pd ~ ztc {no + 2)b'^d + ^(no + 2){no + l)b'^a^ ^ 

since b,d ^ 1 and SE2h — ztc — U{no + 1) — /i > 0. 
Thus near ts — tc, the energy cost from _Eq„ clearly 
outweighs energy gain from E[^ (since ztc/U <C 1 for 
all parameter regime and b,c,d,e ^ 1 near ts — tc) 
and hence the coefficients d and e remain vanishingly 
small at the phase transition. As we move inside 
the superfluid phase and ts becomes large compared 
to tc, these coefficients becomes significant. In the 
present work, we have always restricted ourselves to 
regions sufficiently near the transition line where the 
coefficients d and e are small. We note that we have 
explicitly checked numerically from the full variational 
energy ( including 0{t\/U'^) terms which we have not 
written down explicitly here to avoid clutter), that the 
above-mentioned qualitative argument always holds for 
all the regimes studied in the present work. We shall 
therefore neglect these terms in the rest of this section. 

To compare the results from different approaches, we 
now plot the phase diagrams for ng = 1 obtained by min- 
imizing the full variational energy E^ to 0(t^/L/^)(Eq. 
IA1[) , from the mean-field equation (Eq. IA5[) , and by com- 
puting the energy of the defect states to 0{t'^/U^)^. 
These plots are shown in Fig. [9] for d — 2 and in Fig. 
[TUl for d = 3. We find that in spite of the dissimilar- 
ity of the two approaches discussed above, the numeri- 



cal phase boundary obtained by minimizing Ey matches 
that obtained from 0{t\/U^) defect state calculation of 
Ref. m quite well in both cases, but differs substan- 
tially from the phase boundary obtained using the mean- 
field theory. Further, for d = 2, Quantum Monte Carlo 
data, available for the critical hopping at the tip of the 
Mott \oh(M^, predicts {tc/U)MC = 0.061 ± 0.006. The 
corresponding values obtained from minimization of E.^ 
and 0(tg/?7^) defect state energy calculation of Rcf. 
m are (tc/U)^^, = 0.0755 and (ic/f/)dotoct = 0.0735 
respectively. The mean-field theory predicts a value 
{tc/U)^f = 0.041 while an 0{t%/U^) calculation of de- 
fect states energies^^ gives (tc/C^)dcfoct = 0.068. Thus 
we conclude that the phase boundary obtained from our 
variational energy calculation scheme compares well with 
the defect state energy calculation to 0(t^/J7^). 

In 0? = 3, we can compare our results in Fig.[TU]to Figs. 
[5l [6] and [7] for A = and /i = 0.7 and /i = 0.4 respectively. 
We find that the phase diagram obtained using the pro- 
jection operator techniques predicts tc/U — 0.0232 for 
/i/t/ = 0.7 and = 5 (Fig.© and predicts tc/U = 0.0431 
for ^/U — 0.4 and t] — (Fig. [7|). As seen from 
Fig. [ini these values compares well to those obtained 
from defect state calculations to 0{t\/U^). However, 
they do not compare favorably with t™^ which predicts 
tff/[/ = 0.0208 for ^l/U = 0.7 and V^^^ /U = 0.0266 for 
n/U — 0.4. Further, the corresponding values of tc/U 
from the third order defect state calculation are 0.0225 
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pares quite favorably to the projection operator method, 
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obtained from different methods become most appar- 
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